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Abstract 

A model for verifying a consistency of the next-to-leading order hadronic contributions to 
the muon anomalous magnetic moment with those of the leading order is proposed. A part 
of the next-to-leading order hadronic contributions related to the vacuum polarization is 
rather accurately reproduced in the model. I find a new numerical value for the light-by- 
light hadronic contribution that leads to agreement with recent experimental result for 
the muon anomalous magnetic moment. 



1 Introduction 



A numerical value of the muon anomalous magnetic moment (MAMM) measured experimen- 
tally with high precision can be used to quantitatively test the theories suggested for describing 
particle interactions. The experimental result for MAMM presented in ref. |I[ reads 

a;^^P = 116 592 023(151) x 10"^^ (1) 

with the uncertainty 151 x 10~^^. The main anomalous effect is due to Schwinger term 

= f (2) 

where a is the fine structure constant a"^ = 137.036 . . . The theoretical contributions presently 
computed in the standard model for the comparison with the experimental value given in eq. (j^) 
are divided into three parts: leptonic (QED), electroweak (EW), and hadronic (had) one. The 
pure leptonic part is computed in perturbative QED through order p|, ||. The numerical 
value of the QED contribution to the muon anomalous magnetic moment reads (as a review 
see i) 

aj^™ = 116 584 705.7(2.9) x 10"^^ . (3) 

The EW corrections are well defined in the perturbation theory framework of the standard 
model and have been computed with the two-loop accuracy (as a review see [^) 

a)^^ = 152(4) X 10-" . (4) 

Numerically this contribution matches the present experimental uncertainty. The EW cor- 
rection will be noticeable if a goal to reach the planned experimental accuracy 40 x 10"^^ is 
accomphshed (as a review see ||l[|). 

The hadronic contribution to MAMM is sensitive to the infrared region and cannot be 
computed in perturbative QCD with light quarks. The current masses of light quarks are too 
small to provide a necessary infrared cutoff and explicit models of hadronization are required 
for the quantitative analysis. This constitutes a main difficulty of the theoretical analysis of 
MAMM in the standard model. Writing 

aT = ar + «r + «r (5) 

and assuming 

= «r (6) 

one finds a numerical value for the hadronic contribution to MAMM 



had 



Ith = (7165.3 + 151|exp + 2.9|qed + 4|ew) x 10"" (7) 
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with the experimental error dominating the uncertainty. 

Since the hadronic contribution is sensitive to the details of the strong coupling regime of 
QCD at low energies and cannot be unambiguously computed in perturbation theory frame- 
work the theoretical prediction for MAMM in the standard model depends crucially on how 
this contribution is estimated. In the absence of reliable theoretical tool for computation in this 
region one turns to experimental data on low-energy hadron interactions for extracting a neces- 
sary numerical value. In general terms the hadronic contribution to MAMM is determined by 
the correlation functions of electromagnetic (EM) currents. As the source for the EM current 
is readily available for a wide range of energies one tries to extract these functions (or some 
their characteristics) from experiment. Without an explicit use of QCD the correction is 
generated through the EM interaction ei^°''^A^ with jj"'^ being the hadronic part of the EM 
current. At the leading order (a^ in the formal power-counting) only the two-point correlation 
function of the EM currents emerges in the analysis of hadronic contributions to MAMM 

~ (j;"'^(x)jr (0)) . (8) 

At the next-to-leading order (a^) the four-point correlation function appears 

n4 ~ (jj'''(x)jr (i/)j''^'(^)jr(o)) • (9) 

These correlators are not calculable perturbatively in the region that is essential for the deter- 
mination of the hadronic contributions to MAMM. The leading contribution to MAMM comes 
from the two-point correlator eq. (P) referred to as the hadronic vacuum polarization contribu- 
tion while the four-point function eq. (^ first emerges at the order, most explicitly as the 
light-by-light scattering. To avoid using QCD in the strong coupling mode one has to extract 
the necessary contribution to MAMM by studying these two correlation functions experimen- 
tally without an explicit realization of the hadronic EM current in terms of elementary 
fields. Historically this was a way of studying the EM properties of hadrons before emerging 
QCD fundamental theory of strong interactions (e.g. [H]). 



2 Hadronic contribution at the leading order 

At the leading order in a the hadronic contribution is described by the correlator 

I j {Tf^'^\x)jl:^'iO))e''^^dx = {q,q, - g ,,q^)I]^-''\q^) (10) 

which reduces to a single function I\^^{q'^) of one variable q^. The correlator is transverse 
due to conservation of the hadronic EM current in the standard model. This function gives a 
contribution to the muon anomalous magnetic moment (e.g. p) 

a^fihO) =Ati(-X r —K{s)\m n^^'i(s) (11) 

^ Vtt/ Jiml S 
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with the one-loop kernel of the form 

K{s)= f'dx ,f[!~''| . ■ (12) 
Jo + (1 - x)-^ 

^ ' raj, 

Here Im n'^^'^(s) = Im {n^'^'^(g^) |^2=^_,_jo}, is a muon mass. 

The leading order hadronic contribution to MAMM is represented by an integral over the 
hadron spectrum. No specific information about the function Im n^'^'^(s) is necessary point- 
wise except its threshold structure in the low-energy region. For the applications at the leading 
order in a the function Im 11^'^'^ (s) can be uniquely identified with data from e^e~ annihilation 
into hadrons. Introducing the experimental ratio 

g.P(,) , - hadrons) (13) 

and identifying it with the theoretical quantity B}^\i^o{,s) taken at the leading order in a as 

= 127rlm n^^"'^(s) (14) 



one finds 



"HLOj^j^ (15) 



The contribution to MAMM based on the representation given in eq. (^5|) is well studied. 
Several recent determinations are 

aJ"'^(LO) = 7011(94) x 10"" (ref. [6]); (16) 

aJ^<i(LO) = 6924(62) x lO"" (ref. [7]); (17) 

aJ'^d(LO) = 6988(111) x 10"" (ref. [8]) . (18) 

I use in the further analysis a naive average of these three results (both central values and 
errors are averaged) which reads 

aJ^'^(LO) = 6974(89) x 10"" . (19) 

Writing 

and comparing with eq. (0) one has (in units 10^^^) 

aJ^<^(NLO) = 7165 + 151|exp + 2.9|qed + 4|ew - 6974(89)|lo 

= 191 + 151|exp + 2.9|qed + 4|ew + 89|had • (20) 
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Assuming the statistical independence of the uncertainties one finds after adding them in 
quadratures 

aJ^'^(NLO) = (191 ± 175) x 10~^^ (21) 

that does not allow to see higher order hadronic effects clearly. The error comes mainly from 
the experimental value eq. (p and the leading order hadronic data eq. (^) the statistical 
correlation of which is supposed to be small as they come from different sources. Other errors 
are negligible. For the target experimental error of MAMM at the level of 40 x 10"^^ one finds 
the uncertainty of the NLO hadronic contribution to become 98 x 10^^^. Assuming that the 
mean value of a^^^ in the planned experiment will not change one finds a numerical value for 
the NLO hadronic contribution 

aJ'^'^(NLO) = (191 ± 98) x 10"^^ (22) 

that makes the NLO hadronic effects noticeable at the level of two standard deviations. If the 
mean value of a^^^ will change in the range of the present experimental uncertainty 151 x 10~^^ 
the NLO hadronic effects can be more or less pronounced. From the naive counting in a 
a numerical value for the theoretical NLO hadronic contribution about 50 x 10~^^ can be 



expected. This number is comparable in magnitude with the uncertainty in eq. (22) and should 
be taken into account. 



3 Hadronic contribution at next-to-leading order 

In NLO there is no such a transparency with determining hadronic contributions as in LO. 
Basically there are two new features. On the experiment side the interpretation of data to be 
used in the NLO theoretical calculations is more involved. The problem is to avoid double 
counting as a part of the hadronic contributions has already been accounted through the use of 
data at LO. On the theory side a new correlation function 114 from eq. (^ which is much more 
complicated than the two-point correlator enters the game. At present there is no accurate 
experimental determination of the four-point function in the kinematical range necessary for 
MAMM computation and one has to rely on phenomenological models used for this function. It 
is difficult to control the accuracy of such models that introduces an explicit model dependence 
in the calculation of the NLO hadronic contribution and makes predictions less definite than 
in LO. 

3.1 Interpretation of data at NLO of EM interaction 

For applications at NLO in a the extraction of data is more involved. For instance, one 
should explicitly take into account the NLO corrections to theoretical factors that emerge in 
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a description of the process from which a set of data is taken. These "theoretically corrected" 
data should be used in NLO applications for computing MAMM. As the sets of data are mainly 
extracted from e~^e~ annihilation I discuss this particular process in some detail. 

3.1.1 One-photon mediated e+e annihilation 

The main object studied experimentally in this sector is the full photon propagator D{q^) 

with n(g^) = n''^P*(g^) + Yi^'"^{q^) being a one-particle irreducible block, = Ana. Note that 
in higher orders the one-particle irreducible block does not split into a sum of pure leptonic 
and pure hadronic contributions. It happens first at NNLO which is far beyond the practical 
interest though. I discuss only NLO or terms in the formal a power-counting. Since the 
data are collected at low energies the EW sector can be excluded. With these restrictions the 
cross section of e^e~ annihilation through the one-photon exchange at NLO without vertex 
corrections to initial states is proportional to 

e^Imllfs) 



Im {g D{q )|,..,^.o} - (i + e2Ren(.))2 + e^Imn(s)^ 

e2(Imn'*^P*(s) + Imn^^'i(s)) ^^4) 



(1 + e2Ren(s))2 + e^\mli{sf 
The theoretical expression for the R ratio at NLO reads 



, _ imn^-'^(s) 
^ '^^^^'^ - imnM.(s) ■ ^^^^ 



If R \nlo{.s) is identified with R^^^[s) from eq. ([TB|) then Imll (s) can be restored by using 



a theoretically calculated Imn^'^(s). For s ^ one finds with NLO accuracy 

127rlmn^"'^(s) = i?'="P(s) (l + • (26) 

In some analyses the cross section a{e^e~ —>■ hadrons) divided by the normalization factor 

= ^ (27) 

is used as a data set Then the relation 

127rlmn'''''^(s) = a(e+e- ^ hadrons) /(Tq + 0(a) (28) 

is valid only at the leading order in a. One of the differences with the R ratio at NLO from 
eq. (|25| ) is the term Ren'^'^'^(s) from the denominator in eq. (p^). The quantity Ren'^^'^(s) 
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can be found by reiterating the leading order term Imn'^'^'^(s) through the dispersion relation 
that gives a relative error of o? order. The NLO contribution in the denominator is related 
to the running of the EM coupling constant and can partly be taken into account through the 
renormalization group technique for the energies far from the resonances 0. Another difference 
is the corrections to the production vertex that should properly be taken into account as they 
enter the cross section. Extracting Imn'^'^'^(s) from the cross section requires to subtract these 
corrections from the data in order to avoid double counting in the NLO analysis of MAMM if 
a theoretical NLO kernel for averaging the two-point correlator is used. 

The use of the i?-ratio is preferable from the theoretical point of view as it relates data 
to the imaginary part of the two-point hadronic correlator Imn'^'^'^(s) in a simple way. It is 
also preferable from the experimental point of view since the total normalization of the data is 
fixed that helps to eliminate systematic errors. In this respect a r-data set that can be used 
to determine that part of the two-point function that is generated by the isovector part of the 
hadronic EM current in the limit of exact isotopic invariance has a different normalization at 
NLO and should be corrected by an explicit account for contributions of the relative a order. 
Note that the NLO correction emerging from the interpretation of data can be controlled 
theoretically within counting in a while corrections due to the violation of isotopic invariance 
for the data obtained from the r and e^e~ channels can only be estimated in models. The 
problem of different normalization remains also for heavy hadrons if their contribution to the 
cross section is calculated from their leptonic branchings. The NLO contribution of heavy 
flavors is not essential though because it is small. 

3.1.2 Two-photon mediated e^e~ annihilation 

The NLO cross section of e^e~ annihilation contains a contribution of the two-photon anni- 
hilation with one hadronic insertion into the photon propagator. This contribution requires 
a special treatment before the data set is related to the hadronic two-point function with the 
NLO accuracy. For instance, the NLO kernel for the MAMM diagram with a vertex correction 
integrates the part of data emerging through the double photon scattering channel in e^e~ 
annihilation. This can lead to double counting at NLO for MAMM. 

Thus, one concludes that at NLO a strict correlation between sets of data and theoretical 
expressions for the NLO kernels emerges. This leads to additional contributions of the relative 
a. order that numerically amounts to about 1% of the leading order contribution which is the 
precision one is trying to reach for the comparison with the experimental result for MAMM. 
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3.2 Four-point correlator 

At NLO a new correlation function of hadronic EM currents emerges. This correlation function 
leads to a new effect which is known as light-by- light scattering. Besides this explicit effect a 
less pronounced mixed effect also emerges. The four-point function gives a contribution to the 
full photon propagator of the form 



where D^^{x) is a free photon propagator with a scalar amplitude D{x) ~ In other words 

a projection of the four-point function of the form 



is present in the two-photon Green function. In QCD and other models where the EM current 
is explicitly expressed through the elementary fields this contribution is interpreted as an EM 
correction to the one-particle irreducible block. There is an option to include this contribution 
to the two-point function. I do not consider this option since the picture of the local interaction 
of the photon with the hadronic EM current is lost in such case. 

Thus, an accurate account of NLO hadronic contributions to MAMM from general principles 
is rather a challenging task both experimentally and theoretically. As a first approach to 
it one can use an effective theory with few free parameters providing a unique framework 
for calculations at LO and NLO. In such an approach the LO information is used to obtain 
numerical values for the model parameters. The NLO results are then computed theoretically. 
This approach can also serve as a base for verifying a consistence of the estimates for the NLO 
hadronic contributions made in different phenomenological models. 

4 A model for hadronic contributions 

In this section I describe a model to check a consistence of the NLO hadronic contributions 
and especially the light- by-light contribution with the results of LO analysis for MAMM. The 
simplest version of the model contains three light quarks with QCD quantum numbers and the 
mass nig which is the only model parameter. The numerical value of is fixed from the LO 
hadronic contribution and then used to find the NLO result. Heavy quarks enter the model 
with their standard masses. In this model the calculations are explicit and can be performed 
analytically that is an advantage. Indeed, the model differs from the leptonic sector only by 
the QCD group factors and the numerical values of fermion masses. 



(29) 




(30) 
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4.1 Fixing ruq from the LO hadronic contribution 

A fermion with mass rUq without QCD group factors (as a lepton) gives the LO contribution 
to MAMM of the form 

af-(LO) = /K)(^)' (31) 

with 

= r PMlMds (32) 

and 

Exphcit integration over s with the kernel K{s) from eq. (|T^ gives 

I{mq) = / dx{l — x)[—'K{x,mq)] (34) 
Jo 

where 



and 



7r(x, rriq) = Q- - 1^ y^{z) - ^ (35) 

<^(z) = ^ArcTanh(v^) - 1, z= -—- ^ — . (36) 

y'z 4m^(l — x) + m^x^ 

An analytical expression for the function I{mq) is known, however, the integral representation 
given in eq. (p^ ) is sufficient for practical applications. 

Contributions of the heavy c and b quarks can directly be computed in QCD perturbation 
theory independently of the model. In the present calculation I use only free quark approxima- 
tion for simplicity. For the c quark with the mass rric = 1.6 GeV and the charge Cc = 2/3 
one finds from eq. (|3lD multiplied by the group factor 3e^ = 4/3 

07^^(10; c) = 69.3 X 10-". (37) 

The 6-quark contribution for rrif, = 4.8 GeV [|1T| and ej, = — 1/3 is small and reads 

a7<^(L0;b) = 1.9 X 10^". (38) 

Thus the contribution of light hadronic modes that is represented in our model by light fermions 
with the mass amounts to 

a;;;^°'^(LO; uds) = (6974.3 - 69.3 - 1.9) x 10"" = 6903(89) x 10"" . (39) 

I assume that this result directly corresponds to the contribution of the two-point correlator 
at the leading order as given in eq. (1TT|). It means that a real data set is properly corrected to 



9 



extract rmlmll^^'^{s). As was discussed above the extraction of r'mImU^''^'^{s) with the NLO 
accuracy requires a careful interpretation of data which is assumed to be done. 

From eqs. ( pl| - |36|) and eq. (|39|) one obtains a numerical value for the single model parameter 



rriq 

rrig = 179 ± 1 MeV . (40) 

This value is rather close to that of the charged pion mass, = 139.6 MeV, which is expected 
as the LO contribution is mainly sensitive to the first derivative at = of the two-point 
function n'^^'^(g^) given in eq. (0). For definiteness I give the LO contribution of the light 
hadronic modes within the model obtained literally with the number from eq. (^0]) 

a;^°'^(LO; uds)(mg = 179 ± 1 MeV) = (6920 ± 70) x 10"" . (41) 

Since in the framework of the model the NLO hadronic corrections to MAMM are determined 
by the single parameter rriq with the numerical value from eq. (^0[) they can readily be found. 

4.2 Hadronic contributions at NLO 

The first check is to use the model for computing the higher order hadronic corrections due to 



vacuum polarization graphs. The data-based analysis gives for the NLO effects of this type [12 



aj^'^(vac; NLO) = -101(6) x 10"^^ . (42) 

This number is about 1.5% of the leading term as expected. As was discussed above at this 
level of precision the numerical value for the NLO contribution depends strongly on the data 
sets used in the analysis. For different data sets the different expressions of the NLO kernel 
should be used to avoid double counting. For example, if the i?-ratio is used in the one-loop 
computation then the LO result should first be divided by the factor (cf. eq. (^61) ) 

127rlmn'^^(.)|,»„. = {^ + 1^ + (^)) (43) 

before being used in the NLO analysis that changes the LO result by 12 x 10^^^ that exceeds 
the uncertainty quoted in eq. (^2]). In fact, even mass suppressed terms can be important at 
this level of precision and the entire function Imn'^^(s) should be integrated since the mass 
terms from the leading order can partly cancel the NLO corrections in a. For other types of 
data (r data especially) the change can be larger. This uncertainty is a reflection of the mixture 
of contributions at NLO. 

In the proposed model the analysis is unambiguous and straightforward. I present different 



contributions separately for a detailed comparison with the results from ref. |12 |. 

For the vertex type contributions I use the explicit analytical formulae in the leading order 
of the mass expansion as they are given in ref. W^. The exact expressions are presented in 
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ref. ||T^ . The analytical expression for the contribution of a fermion with mass ruq without any 
group factors reads 



fprm/ \ 2/m„\^/ 2689 vr^ 23 , m„ , , ^ , 

term/ m„) = — — \ \ In — ^ - (44) 

^ ' ?>\mj \ 5400 15 90 mj\n} ^ ' 



a 
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that leads to a numerical result for the light mode contribution in the model 

a;;^°<^(ver; NLO; uds) = -172 x 10"^^ . (45) 

A more accurate evaluation (using numerical integration with the kernel given up to the third 
order in the mass expansion from ref. [Q) gives for the contribution of light modes 



^mod^^g^. ^LO; uds) = -188 X 10-1^ . (46) 
The difference with the result obtained by using only the first term of mass expansion given 



in eq. (P5| ) is of order 10%. It is smaller than one could expect from the numerical value of 
the expansion parameter {m^/mqY = (0.106/0.179)^ = 0.36. For the c-quark contribution one 
finds 

a)^°'^(ver; NLO; c) = -4 x 10"" (47) 
while the fe-quark contribution is small 

a7'^(ver; NLO; b) = -0.2 x 10"" . (48) 

The total vertex contribution computed in the model 

a7'^(ver; NLO) = -192 x 10"" (49) 



should be compared with the result of the data-based analysis from ref. Wl 



f- 1^2] (ver; NLO) = -211(5) x 10"" . (50) 



Next check of the model is done for a mixed contribution of the lepton-hadron type. This 
contribution contains the electron and r-lepton loops and depends on three masses m^, rrig, 
and rrif. or m^. For fermions without group factors the contribution is given by the integral 
representation 

a|,^™(db; fi&f2) = ^— ^ j (ix(l - x)7r(x, m/j7r(x, m/J . (51) 
For the combined contribution of light modes with the electron loop one has 

a7'^(db; NLO; e&uds) = 105 x 10"^^ (52) 

and with the r-lepton loop 

a7<^(db; NLO; r&uds) = 0.05 x 10"^^ . (53) 
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The contribution of heavy modes is only visible for the combined insertion of the c-quark loop 
and the electron loop 

a7'^(db; NLO; e&c) = 1.1 x 10"" . (54) 



The results given in eqs. ( |52| - p^ are in good agreement with the data-based estimate 

af- [^2](iept&had; NLO) = 107(2) x 10"" . (55) 

Next comes the contribution from the reiteration of hadronic insertions. The light modes give 

a;^°^(db; NLO; uds&uds) = 3 x 10"" . (56) 

The combination of the c-quark insertion with the light mode loops gives 

a7^(db; NLO; c&uds) = 0.1 x 10"" (57) 

while the contribution of the two c-quark insertions is negligible. The results of the model from 
eqs. (1^,^) are in agreement with the data-based estimates 



a'^^- I^2](had&had; NLO) = 2.7 x 10~" . (5^ 



Thus one sees a good agreement of the model results with calculations based on data. However, 
in the model there is a contribution which is missing in the explicit calculations based on data 
as it is related to the internal structure of the hadronic block. In the data-based calculation this 
contribution is hidden in data while in the model it can explicitly be resolved as a correction to 
the one-particle irreducible hadronic block. At the leading order of the mass ratio the analytical 
expression for this contribution without group factors reads 

<'"(^;NLO;™«)^^(!^)^(^)'. (59, 

The result for the light modes of the model is 

a7<^(4; NLO; uds) = 25 x 10"^^ (60) 
while the c-quark contribution is small 

07*^(4; NLO; c) =0.3 X 10"^^ (61) 
The result for the total NLO hadronic contribution of the vacuum polarization type is 

a7^(vac; NLO) = -58 x 10"^^ . (62) 

The difference with eq. (^2]) comes mainly from two sources: the vertex type contributions 
and the new term related to the one-particle irreducible hadronic block. Both contributions are 
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of the {m^/mqY order that explains the magnitude of the difference. All remarks about the 
double counting in the data-based approach apply here. Only first few terms of expansions in 
the mass ratio {m^/niqY were used for numerical estimates that provided a sufficient accuracy. 

Thus, the model reproduces rather accurately the results for the NLO hadronic contributions 
found in the data-based analysis for the graphs related to vacuum polarization. This has been 
expected as these results are obtained by the integration of the two-point function with the 
NLO kernel. 

The next try for the model is the computation of the light-by-light contribution which is 
given by the four-point correlator. The analytical expression for a contribution of the fermion 
without group factors through the {m^/mqY order reads 



«r,..;NLO;,„,)^g)'{g)^gc(3)-15) 

+ h<\ _ mm 13^(3, - i5L,= - ; . (63) 

\mj \ 810 \mf,J 48600 \m^J 18^^^ 9720 972000y ' ^ ^ 
With this formula one finds the value for the light modes 

a^'^Obl; NLO; uds) = 140.5 x 10"^^ (64) 

and one for the c quark 

a;;^°^(lbl; NLO; c) = 2 x 10"^^ . (65) 
The total light-by-light contribution predicted by the model 

a';^°'^{M; NLO) = 140.5 + 2 = 143 x 10"^^ (66) 

is different from the number used in the literature |jl| 

aj^'^(lbl; standard) = -85(25) x 10^^^ . (67) 

I postpone a discussion of this point till sect. (H). 

Thus, the NLO hadronic contribution obtained in the model reads 

a7'^(NL0) = (-58 + 143) x 10"^^ = 85 x 10"^^ . (68) 

It agrees with the present experimental result from eq. (|2l| ) which we repeat here 

al""^ = (191 ± 175) X 10"" . 

The agreement with the future experimental result for MAMM depends on a possible change 
of the mean value of a^^^ as one sees from eq. (p^). 
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The prediction of the NLO hadronic contribution obtained in the model is fairly sensitive to 
the numerical value of the mass parameter for the light modes. This numerical value is however 
strictly determined by the LO result. To check how sensitive to the details of the model the 
results are I introduce a mass difference between s and u, d quarks (S'f/(3) flavor violation in 
the mass sector in the approximation of exact isotopic invariance). I write = iriq + 0.18 GeV 
with 0.18 GeV being the value of the running mass for the strange quark. Then one finds 

TRq = 166 ± 1 GeV (69) 

with 

^modi^LO; uds)(mg = 166 ± 1 MeV) = (6928 ± 71) x 10"" . (70) 

The prediction of the NLO contribution in this case as compared to the SU (3) symmetric one 
is 

^modi(j^LO) - a7'^(NL0) = (2 + 4 + 14) X 10"" (71) 

where the first term comes from vertex corrections, the second term comes from insertions into 
the photon propagator and the last term comes from the light-by-light graphs. The result is 
fairly stable. Finally, the model with SU{'i) flavor violation in the mass sector gives the NLO 
hadronic contribution to MAMM 

^modi^^LO) = 105 X 10"" (72) 

which is rather close to the prediction of the model with S'?7(3) symmetric mass arrangement 
from eq. (p^). 

One could consider an even more sophisticated model including a violation of the isotopic 
invariance by using different masses for u and d quarks. An additional uncertainty emerges from 
the errors in the numerical value for the c-quark mass. By using the MS mass around 1.3 GeV 
for the c-quark one could enhance its LO contribution by about 50% (a leading order rescaling 
factor is (mc(pole)/mc(MS))^ = (1.6/1.3)^ = 1.5). Within the proposed model the use of the 
pole mass of the heavy quark looks more natural while an account of the difference between the 
numerical values for the pole and MS masses is beyond the accuracy of the approximation used 
for heavy quarks. It can readily be done since the contribution of heavy quarks is perturbative 
and corrections in the strong coupling constant can reliably be found. 

5 Discussion 

The underlying idea of the presented analysis is to introduce a framework for computing the 
NLO hadronic contributions to MAMM using the LO information. Presently the results for the 
light-by-light contribution that is the most interesting term at NLO are available analytically for 
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fermions that dictates the choice of the model from the technical point of view almost uniquely. 
Thus, a model of massive quarks with the EM interaction emerges as a suitable candidate. It 
is not an approximation for QCD as a gauge model with constituent quarks. It is just a bridge 
from LO to NLO results for hadronic contributions to a particular observable. Note that for 
another important parameter of the standard model - the running EM coupling constant at 
the scale of the Z-boson mass - there is no possibility to use such kind of a model as there is no 
important next-to-leading order terms to compute. Calculations for the infrared (IR) sensitive 
observables using constituent quarks with masses around 300-500 MeV as the only IR scales 
are unjustified in pQCD in general since the higher order corrections in the strong coupling 
constant cannot be found. Also the introduction of finite masses for the light quarks explicitly 
violates chiral invariance which is a well established symmetry of the light hadronic sector. 
In this sense the approximation for QCD with constituent quarks cannot be considered as a 
reasonable general framework. In the high energy limit the massless approximation for strong 
interactions is perturbative and quite precise. This means that high energy contributions to 
MAMM can be represented by almost any model that satisfies the duality constraints. In this 
sense the fermionic model fits the standard approximation for large energies. However, the main 
contribution to MAMM comes from the IR region where there is no sensible approximation for 
strong interactions deduced from QCD. Therefore the necessary characteristics of the strong 
interaction amplitudes relevant for the computation of MAMM has to be extracted from data. 
The first amplitude that emerges is the two-point correlator that is given by a single function 
of one variable with simple analytic properties (see eqs. (|8| JT0|) ) . For computing MAMM one 
need not know the point-wise behavior of the spectrum but only the integral over all energies 
with some enhancement of the threshold region. A model of massive fermions is then well 
suitable to fit this integral over data. When hadrons are introduced into the threshold IR 
region to fit experiment the effective masses of quarks increase. Therefore an account of low- 
energy hadronization for the two-point function entering MAMM is achieved by introducing an 
explicit cut in energy in the sum over the states. In practice, at the leading order the hadron 
contributions are represented by the pion with an EM interaction of the form ej'^A^ where at 
the leading order = i{n^ d^i tt^). The inclusion of pions leads to the scalar type of the 
spectrum near the threshold 



instead of the fermionic form given in eq. (|33D. Furthermore, the fermionic contributions can 
be moved to higher energies by using vector mesons. In the vector meson dominance model 
one identifies the EM current with the canonically normalized elementary p-meson field 
through the relation = fpP^. Here fp gives a form factor related to the leptonic width of 
the p-meson. Because of the nature of the MAMM observable this contribution can well be 




(73) 
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represented by fermions already as it resides at a rather large scale. This hadronization picture 
is transparent for the two-point function which is sufficient for the LO analysis. At NLO an 
hadronization procedure for the four-point function is necessary. Within a hadron picture of 
the low-energy spectrum the most important contribution to MAMM comes from pions. To 
quantitatively handle contributions from the four-point function a quantum field model for 
pions given by the Lagrangian 

Llowenergy = \D^^Tlf - m^TT^ D^ = d^- 16 A ^ (74) 

is introduced. This model generates vertices that allow to compute the pion contribution to 
the four-point hadronic EM current correlator that enters the light-by-light diagram explicitly. 
The high energy contribution of this model should then be replaced by the standard quark 
contributions. In the pure fermionic model with the small effective mass this replacement is 
effectively made at rather low energies that makes the separate contribution of pions small 
or even vanishing. Thus, the hadronization procedure of the model is realized through light 
massive quarks rather than real hadrons. Note that the hadronization picture need not be 
universal for all strong interaction processes but can specially be tailored for a given observable. 

The results for MAMM related to the two-point function which have been obtained in data- 
based analysis arc well reproduced by the model with the mass of the light fermion around the 
pion mass. Using the model prediction for the light-by-light graph I find an agreement of the 
NLO hadron contribution to MAMM with experiment. The results for the light-by- light graph 
in the pion model arc available numerically. In the absence of analytical expressions for the 
light-by-light contributions in the pion model I could not quantitatively check how fermionic 
contributions replace the pion ones when the effective fermion mass decreases. However, it 
seems probable that the explicit inclusion of the pion contributions in the framework of the 
present model will shift the effective mass of light quarks larger. 

The fermionic model gives a smooth spectrum at low energies. An important question is 
whether such a smooth spectrum is a reasonable approximation for computing MAMM. In the 
two-point correlator there are no resonances in the relevant region. The contribution of the 
p-meson is located at relatively large scales. In the axial channel, for instance, the situation is 
different because of the presence of the pion resonance and probably the model with massive 
quarks would not fit. Note, however, that as soon as the pion is considered to be massive (not 
a pure Goldstone mode) the chiral invariance is explicitly broken that makes quarks massive 
as well (or vice versa). Note also that a model can be suited for a description of a specific 
observable and need not give a universal approximation of any Green function. For instance, in 
the axial-vector two-point correlator the projection related to spin-one particles contains only 
massive resonances and the spectrum can well be approximated by a fermionic model without 
the pion pole. For the four-point functions the situation is more complicated though. In the 
literature there are models where the four-point function at low energies is represented through 
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the elementary fields of neutral pseudoscalar bosons in order to compute contributions of the 
light-by-light graphs to MAMM. The representation employs the neutral pion contribution to 
the four-point function through the iteration of an effective Lagrangian for the interaction of 
the neutral pion 7r° with photons due to the Abelian anomaly in axial current (as a review 



see ref. ||T6[). The result for the light-by- light contribution obtained using the neutral pion 
dominance eq. (^) is different from one obtained in the present model eq. (^). In general, the 
reduction of the four-point amplitude of the hadronic EM currents to a two-point correlator of 
axial currents uses the operator product expansion at small distances 

^2^J^'(^)jr(-^)U^o = e,,^,x-j^{0)C{x') + ... (75) 

where C(x^) is a coefficient function of the local operator ^5(0) which has quantum numbers 
of axial current (e.g. [l^). In other words the combination of two hadronic EM currents of the 
form 

e,u.xeTf,^'{x + Oj-r (x - OF{e) (76) 

taken at small ^ with some form factor -F(^^) may act in some applications as a local axial 
current that can serve as an interpolation field for the pion. Thus, this combination can be 
replaced by a fundamental pion field in a hadronization procedure. This kind of factorization 
for the four-point amplitude is valid for 77 — > 77 scattering in a specific region of the phase 
space in kinematic variables where all three external momenta are essential. In other regions 
of the phase space the saturation of the scattering amplitude with the pion pole contribution 
is invalid. The projection of the four-point function that emerges in the light-by-light graphs 
for the MAMM calculation has the form 

dxx^(Tj;'^'^(x)jr(y)jr(^)jJ'^'^(0)) . (77) 

In momentum space this projection depends on two external momenta only as the third mo- 
mentum is set to zero after differentiation according to the definition of MAMM. In the neutral 
pseudoscalar model the projection of the four-point function given in eq. ([77|) is saturated by 
the contribution of the neutral pion that seems to be invalid in the kinematical region relevant 
for MAMM computation. In the absence of the neutral pion pole contribution in the hadroniza- 
tion picture for the light-by- light graph the fermionic model can be used for its computation 
on the same footing as it was used for vacuum polarization graphs. In fact, the neutral pion 
contribution gives the major difference with the present analysis based on the fermionic model. 
However, the corresponding contribution of the neutral pion to the projection of the four-point 
function emerging in the photon propagator is usually not considered. In other words, the 
neutral pion approximation for the four-point function should also be taken into account in 
eq. (|30D as it is accounted for in eq. ([77|) . If the neutral pion contribution to eq. (|30D does 



not vanish by some symmetry considerations it can lead to a cut starting from the pion mass 
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square that seems to contradict the threshold behavior of the spectrum known in e^e" 
annihilation. This calls for a necessity to evaluate the validity of the neutral pion dominance 
model for the four-point function in the kinematical region relevant for computing the NLO 
hadronic contribution to MAMM. 

Despite the fact that the fermionic model with the mass = 179 MeV predicts the value 
for the NLO hadronic contribution to MAMM in agreement with experiment there remains 
a disturbing feeling that this prediction is obtained within an unrealistic approximation for 
strong interactions and, therefore, cannot be taken seriously. A historic reminiscence may be 
appropriate here. A century ago thinking about the light as existing in the form of discrete 
portions - photon quanta - was rather disturbing for classical physics. However, the quantum 
representation allowed for the quantitative explanation of experimental facts on photoeffect 
and black body radiation. It did not change the description of electromagnetic phenomena 
insensitive to the quantum nature of the light. It may happen that the muon anomalous 
magnetic moment is sensitive to the contribution of all hadrons in a way it would be sensitive 
to that of free fermions with an appropriate mass which is a standard realization of duality 
concept. The direct application of this concept to a particular case of MAMM looks suspicious 
because the IR region is explicitly involved in the analysis and the results depend strongly on 
the numerical value of the effective quark mass which happens to be rather small. The model, 
however, is only designed for computing the NLO hadronic contribution to MAMM using the 
LO result as input. This does not mean that this model approximation suited for computing 
MAMM is in any sense a universal limit of QCD automatically applicable to other observables. 

6 Conclusion 

A model for describing the NLO hadronic contributions to the muon anomalous magnetic 
moment is proposed. The model contains a single parameter that is fixed from the experimental 
result for the LO hadronic contribution to MAMM. The model describes the NLO hadronic 
contributions of the vacuum polarization type in agreement with existing estimates. However, it 
predicts a numerical value for the light-by-light contribution which is different from one used in 
the literature that considerably changes the prediction of the total NLO hadronic contribution 
to MAMM. The prediction of the model agrees with the present experimental value for MAMM. 
A resolution of the contradiction between the estimates for the NLO hadronic contribution, or 
rather for the light-by-light contribution, obtained in the present model and existing in the 
literature could help in verifying the validity of the standard model. 
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